AFM Dissipation Topography of Soliton Superstructures in Adsorbed Overlayers 
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In the atomic force microscope, the nanoscale force topography of even complex surface superstruc- 
tures is extracted by the changing vibration frequency of a scanning tip. An alternative dissipation 
topography with similar or even better contrast has been demonstrated recently by mapping the 
(x, y)-dependent tip damping: but the detailed damping mechanism is still unknown. Here we iden- 
tify two different tip dissipation mechanisms: local mechanical softness, and hysteresis. Motivated 
by recent data, we describe both of them in a one-dimensional model of Moire superstructures of 
incommensurate overlayers. Local softness at "soliton" defects yields a dissipation contrast that can 
be much larger than the corresponding density or corrugation contrast. At realistically low vibration 
frequencies, however, a much stronger and more effective dissipation is caused by the tip-induced 
nonlinear jumping of the soliton, naturally developing bi-stability and hysteresis. Signatures of this 
mechanism are proposed for experimental identification. 

PACS numbers: 46.55.+d, 07.79. Lh, 07.79.Sp, 81.40.Pq, 62.20.Qp 



INTRODUCTION 



The tip-based scanning force microscopes of the atomic 
force microscope (AFM) family constitute perhaps the 
single most important tool bag in nanotechnology. The 
substrate topography is extracted from a map of the oscil- 
lation frequency of a tip, hovering a short distance above 
the surface. Besides the frequency shift however, the tip 
also develops a damping, reflecting a position dependent 
mechanical dissipation. Maier et al. [H, [2] showed re- 
cently that AFM dissipation - whose general occurrence 
has been widely discussed by several groups a decade ago 
but whose potential importance was still underesti- 
mated - is able to map exquisitely delicate features such 
as the Moire superstructure pattern formed by misfit dis- 
locations ( "solitons" ) of incommensurate KBr adsorbate 
islands (Fig. on NaCl(lOO). Surprisingly, the exper- 
imental dissipation map, Fig. [Tp, showed similar or bet- 
ter contrast than the corresponding topographic map, 
with a characteristic reversed contrast (higher dissipa- 
tion at the soliton, where topographic height is minimal 
Q ) . Given also the great importance of achieving newer 
routes toward high-quality imaging, this is more than 
a mere curiosity, and deserves a proper understanding. 
Existing linear-response theory and other approaches to 
AFM dissipation [a, 0] and to general frictional dissipa- 
tion [I3-[l2| suggest a larger tip damping above softer 
substrates, and that provides an initial and valuable clue. 
Local tip dissipation can effectively reveal the underlying 
superstructure, since the local mechanical compliance is 
higher for example at surface soliton lines, where atoms 
sit at metastable positions. However at the relatively low 
AFM oscillation frequencies, the current understanding 



rules out linear response as the chief dissipation mecha- 
nism. A typical energy dissipation as large as 0.01 — 1 eV 
per oscillation can only be accounted for by a hysteretic 
response of the interacting tip-substrate system, as was 
understood by theoretical analysis [l3l - fl5j ]. and demon- 
strated experimentally [3, [13] • Such nonlinear effects 
of hysteresis are most likely involved in the surprisingly 
large AFM dissipation contrast of adsorbate superstruc- 
tures too. Yet, it is unclear how inert systems (such as for 
example alkali halide overlayers) could give rise to hys- 
teretic phenomena, and in particular how they would be 
connected to the presence of misfit superstructures ("soli- 
tons" ) . This is the issue which we address here by means 
of dynamical simulations of the simplest one-dimensional 
model. Our main result is the identification of an unex- 
pected soliton-related hysteretic mechanism. During the 
first part of its swing, the tip can locally drag or push an 
underlying defect - here a soliton portion - causing it to 
jump across a (Peierls-Nabarro) energy barrier. During 
the return journey, the defect follows only sluggishly, and 
remains trapped somewhat longer on the wrong side of 
the energy barrier, thus opening a hysteresis loop. The 
area enclosed in the hysteretic force-displacement dia- 
gram represents a large tip energy dissipation, one that 
can survive down to realistically low AFM vibration fre- 
quencies, a regime where the linear-response dissipation 
is quantitatively irrelevant. This mechanism is likely to 
play a significant role every time a "softness pattern" is 
present, and should be easier to observe for horizontal 
than for vertical tip oscillations. 
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FIG. 1: (Color online) AFM Moire superstructures of in- 
commensurate KBr bilayers islands/NaCl(100), adapted from 
Ref. @]: (a) topography; (b) dissipation. Note the opposite 
phase: dissipation is largest at the soliton, where topographic 
height is minimal 8J. (c) ID simulation model with a rigid 
substrate potential of period a S ut>, a harmonic chain of rest 
length ao (the adsorbate overlayer), and the localized tip po- 
tential ut (here attractive) oscillating horizontally (in exper- 
iments the oscillation is usually vertical). 



THE MODEL 

To emphasize the basic and general aspects of the phe- 
nomena, rather than a realistic model targeted on the 
chemical detail of a specific tip/sample configuration [l3l - 



15l . Il8| , we use the simplest possible model - a tip po- 



tential oscillating over a one-dimensional harmonic chain 
(the overlayer) moving in a rigid incommensurate peri- 
odic potential (the substrate), Fig.QF,. The Hamiltonian 
of the mobile overlayer atoms is 

H = E k + U at ^ t + U sub + U T (t), (1) 

where Ek = y J2i x i 1S tne kinetic energy, 



K 



t/at-at = -r- ^(xj+i - Xi - a ) 2 



(2) 



is the mutual (harmonic) interaction potential, and 
U su b — J2i v i x i) is the substrate potential, which we take 
of a pure cosine form: 



v{x) 



-^sub ^sub 
47T 



cos(fc sub x) 



(3) 



Here ao is the mean spacing between adatoms, K is their 
mutual spring constant, fc su b = 27r/a su b, and a su b is the 



period of the substrate potential 19|, 2(j. The two peri- 
odicities ao and a su b define the coverage 9 = a su b/ao- For 
convenience we take a su b as the unit length, F sub as the 
force unit, and the mass m of the particles as the mass 
unit. (To get a feeling for quantities, the frequency units 

^sub "t- -1 ^ 2 a 7ub 2 should be typical of an atomic vibra- 
tion or a Debye frequency wd, typically 1 THz or more). 
The general lack of commensuration between adsorbate 
and substrate periodicities gives rise to two-dimensional 
misfit dislocations, sometimes called solitons, which form 
a regular superstructure with the beat periodicity be- 
tween the two. Fixed boundary conditions (BCs) are 
chosen in order to prevent the advancing tip to drag the 
entire pattern along, that would occur if, e.g., periodic 
BCs were used instead. 

£/t(£) = J2i u T( x i>t) i s the time-dependent oscillating 
potential describing the tip action on the overlayer. We 
represent the AFM tip as a Gaussian-shaped oscillating 
potential, with u^(x,t) ~ u(x — XT(t)), xt(£) = xt + 
Axcos(wTi), and 



u(x) = At exp(— x 2 /a T ) 



(4) 



Here At represents the repulsive (contact AFM, At > 
0) or attractive (noncontact AFM, At < 0) tip-atom 
interaction strength, <tt is the tip width, Ax and lot 
are the tip oscillation amplitude and angular frequency 
around its central position xt- 

The equation of motion for the i-th overlayer atom is 

mxi = -v'(xi) + K(x i+ i - 2xi) + / T (ii,t) - j&i , 

(5) 

where v'(x) = iF su b sin(fc su bx), and the tip force 



h(x,t) 



d_ 

dx 



ur(x,t) 



(6) 



is given by a straightforward analytical expression. A 
damping force term — is introduced to represent all 
dissipation phenomena which remove energy and allow 
the attainment of a stationary frictional state. 

We integrate the equations of motion ([S]) by means of 
a standard adaptive fourth-order Runge-Kutta routine 
2l| starting each simulation from a stationary fully re- 
laxed overlayer, as obtained by a preliminary relaxation 
of equally-spaced adatoms x^(0) = i ■ ao and x;(0) = 0. 
All simulations are carried out at a nearly commensurate 
coverage 8 = 1.06 = ||, realized by means of a chain of 
N = 107 particles in a region of length L = 100 a su b- A fi- 
nite temperature T is implemented by adding a standard 
Langevin random force to Eq. ([5]), and averaging over a 
long simulation time, usually at least 100 tip-oscillation 
periods. The extreme simplicity of the model allows us to 
extend simulations down to the realistic AFM frequency 
in the MHz range, which requires exceedingly long inte- 
gration times. 
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FIG. 2: (Color online) The linear-response regime, (a) Over- 
layer atom density (dashed), rest positions (circles), in the 
periodic substrate potential (solid), (b) The density of zero- 
frequency overlayer vibrational modes peaks at the soliton 
positions, (c) Dissipated power during tip vibration for a tip 
width <tt = At = a S ut>, weakly attractive potential of ampli- 
tude At = — 10~ 4 , 7 = 0.2 (underdamped regime), K = 5 
(fairly rigid overlayer), and oscillation frequency lot = 10 -4 it 
(roughly 1 GHz). This regime is described well by linear- 
response theory. Note a huge dissipation peak contrast at 
the soliton position relative to the terrace between solitons of 
about 10 4 . 

The instantaneous power drained away by the damping 
term amounts to 

-Pdiss = Y" ±i ■ (jXi) = 7 V x\ = — E c , (7) 

^-^ TO 

* i 

and is thus proportional to the total kinetic energy of the 
overlayer. The power pumped by the tip into the chain 
is 

P T = ^/T(Zi,t)-^. (8) 

i 

While these two quantities fluctuate separately, they 
must of course coincide on average over a period r = 
27t/wt in the dynamical steady state 

P = ~ JdtP T {t) = ~ JdtP diss (t), (9) 

also indicating how the work done by the tip oscillation 
is eventually dissipated entirely by the viscous friction 
term. 

RESULTS: LINEAR RESPONSE AND BEYOND 

Figure [5] displays the dissipation results obtained in 
simulation for a weakly attractive tip potential, very high 



oscillation frequency, and general parameters that fall 
well inside the linear-response regime 0, ■ The lin- 
ear response results show (i) strong dissipation enhance- 
ment at solitons, with P(x) several orders of magnitude 
stronger than in a terrace between two of them, closely 
mirroring the phonon local density of states (LDOS); (ii) 
dissipated power which is proportional to A\, indepen- 
dent of the attractive/repulsive sign of the tip-overlayer 
interaction, i.e. of the noncontact or contact mode of 
the AFM; (iii) absolute dissipation values that are very 
weak everywhere, and dropping with decreasing AFM 
frequency as (wt/wd) 2 - Summing up, the predicted rel- 
ative contrast of the soliton pattern in linear response dis- 
sipation is indeed very large. However, the exceedingly 
low value of realistic AFM frequencies (< MHz) relative 
to microscopic frequencies (~ THz) renders this linear- 
dissipation mechanism entirely academic. 



We reach a more realistic regime by enhancing the 
tip-overlayer interaction strength, while still remaining 
in a moderate-interaction regime representing noncon- 
tact AFM. This new regime is dominated by nonlinear 
effects, where dissipation no longer drops as w^,, but at 
most linearly in wt (apart from logarithmic corrections) . 
Comparison of Fig. [3^ with Fig. [2t shows that, nonlinear 
dissipation is again much larger near the solitons than in 
between them. A two-order of magnitude increase in At 
would in linear regime imply a dissipation increase by a 
factor 10 4 , whereas we find a much larger factor of about 
10 6 already at this large frequency (cjt is here 10~ 4 7r, cor- 
responding to the gigahertz range). Decreasing frequency 
down to realistic AFM values, the increase will become 
gigantic, because the nonlinear dissipation lacks the extra 
power of ujt appearing in the linear-regime dissipation. 
The new element brought in by nonlinearity is mechan- 
ical. A strongly interacting tip is now able to drag, or 
to push, the soliton - a mobile entity - forward or back- 
ward during the oscillation cycle. As the soliton must 
overcome the (Peierls-Nabarro) barrier in order to move, 
its motion is sluggish, and can follow the tip only with 
hysteresis and, as anticipated, hysteresis entails a large 
dissipation. As shown in detail in Fig. [3K,b, the higher 
dissipation point A is found to corresponds to two succes- 
sive Peierls-Nabarro barriers being overcome in the oscil- 
latory process, the smaller dissipation point B to a single 
barrier. The dissipation at point C, where the tip poten- 
tial is unable to "grab" the soliton, is negligible by com- 
parison. The onset of this large-dissipation region, domi- 
nated by hysteresis, is rather sharp. Fig. 2] illustrates this 
point, by showing the average power P at location A, di- 
vided by the linear- regime factor | At | 2 , as a function of 
the tip amplitude |At|- Beyond a value of |At| of order 
4 x 10 -3 , the linear regime behavior is abruptly aban- 
doned, and the dissipation increases rapidly by several 
orders of magnitude. 
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FIG. 3: (Color online) The strong-interaction hysteretic regime, (a) Mean power dissipated in the steady regime by a strongly 
interacting tip (At = —0.01, all other parameters the same as in Fig. [2} scanning the same overlayer (K = 5). (b) Force- 
displacement response at three typical scanning points marked in the left panel. The strong dissipation at the solitons is now 
due to hysteretic jumps of the solitons (kinks) across their Peierls-Nabarro barriers. The dissipated energy in a cycle equals 
the area of the hysteretic loop in the force-displacement plane (shaded). When the potential is strong enough to drag or push 
a soliton (points A, B), this occurs with hysteresis and the dissipation is large. With the selected oscillation amplitude, the 
soliton is dragged across two barriers (point A) or a single barrier (point B), depending on the center of tip position xt- When 
the tip grabs no soliton (point C) there is no hysteresis and dissipation drops. 
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FIG. 4: (Color online) The dissipated power normalized to 
the linear regime scaling, P/|At) 2 , versus the strength of the 
tip-overlayer interaction |-i4t], at point A of Fig.[3]and for the 
same parameters as in Fig. [2] Notice the sharp onset of the 
non-linear regime where the hysteretic mechanisms starts to 
play a role, leading to a strong enhancement of dissipation. 



PREDICTIONS AND DISCUSSION 

Our simulated example strongly suggests that a large 
hysteretic component should be present in the existing 

I — ~~L 

dissipation maps [22] of Moire patterns. More generally, 
hysteretic defect dragging should dominate the AFM dis- 
sipation maps. What are the predicted signatures of this 
mechanism? Our model study suggests two main signa- 
tures. 




FIG. 5: (Color online) The dissipation reduction due to ther- 
mal shrinking of the hysteretic loop, for ljt = 10 -4 it at the 
soliton location A of Fig. [3] (a) The frequency dependence of 
the energy dissipated in one period P2-k/ojt, computed for 
several temperatures, exhibiting a clear reduction of P2n/uiT 
due to a rise in temperature, (b) Detail of the right-side area 
of the force-displacement dependency following two typical 
tip-oscillation periods for T = 10 -5 , solid, and 10 -6 , dotted, 
compared to T — 0, dashed: the hysteretic area is shrink- 
ing due to a randomly anticipated thermally activated barrier 
crossing. 



(i) Abruptness of AFM friction onset with increas- 
ing strength of tip-surface interaction. As suggested by 
Figs. [SHU dragging sets in abruptly only above a cer- 
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tain threshold, which means below a certain tip-surface 
or tip-soliton distance. 

(ii) Anomalously mild (linear with logarithmic correc- 
tions) frequency dependence of AFM friction at finite 
temperature. It is a rather general property of all hys- 
teretic friction phenomena to heal away at sufficiently low 
frequencies, where adiabatic motion allows sufficient time 
to jump thermally over barriers. For instance, thermol- 
ubricity experiments 23-2(| and detailed calculations 
within the Tomlison model [27f | show an average fric- 
tion force F = F c -AT 2 / 3 |log(m;/T)| 2/3 , where v is the 
driving velocity, and A and a are system-dependent di- 
mensional constants. In incommensurate overlayers, the 
soliton nearest to the tip behaves similarly to a Tomlin- 
son particle, as is driven across a Peierls-Nabarro bar- 
rier. In Fig. [5^ we do observe a thermal reduction of 
dissipation, due to a shrinkage of the hysteretic loop, il- 
lustrated in Fig. [5}d. When driving is oscillatory as in 
AFM dissipation, the role of v is taken by At wt • We 
find our data to be compatible with a similar relation 
P(T) =P(0)-A'T 2 / 3 \log{a'u} T /T)\ 2/3 . The parameters 
A' and a' are here related to the effective soliton proper- 
ties (mass, damping, barrier height...), and are nontrivial 
functions of the "bare" model parameters. 

We conclude that AFM dissipation maps of incommen- 
surate overlayer superstructures can in principle achieve 
an extremely high contrast resolution of soliton defects 
relative to commensurate terraces. The most important 
theoretically predicted dissipation mechanism is the non- 
linear dragging or pushing of some local portion of the 
defect, where the large tip damping is associated with 
hysteresis of defect motion. Besides a sharp threshold 
in the tip-surface interaction and oscillation magnitudes, 
this mechanism predicts a very characteristic logarithmic 
dependence (eventually turning to linear at extremely 
low frequencies) of dissipation upon frequency and tem- 
perature. More generally, the nonlinear dragging of soft 
defects or features (e.g. a floppy residue in a biomolecule) 
should give rise to a strong visibility in AFM dissipation 
topography, of considerable potential impact for applica- 
tions. 
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